Membrane transporters in drug development and as determinants of precision medicine.

The effect of membrane transporters on drug disposition, efficacy and safety is now well recognized. Since the initial publication from the International Transporter Consortium, significant progress has been made in understanding the roles and functions of transporters, as well as in the development of tools and models to assess and predict transporter-mediated activity, toxicity and drug-drug interactions (DDIs). Notable advances include an increased understanding of the effects of intrinsic and extrinsic factors on transporter activity, the application of physiologically based pharmacokinetic modelling in predicting transporter-mediated drug disposition, the identification of endogenous biomarkers to assess transporter-mediated DDIs and the determination of the cryogenic electron microscopy structures of SLC and ABC transporters. This article provides an overview of these key developments, highlighting unanswered questions, regulatory considerations and future directions.

Genome-Wide Association Study Identifies Pharmacogenomic Variants Associated With Metformin Glycemic Response in African American Patients With Type 2 Diabetes.

OBJECTIVE: Metformin is the most common treatment for type 2 diabetes (T2D). However, there have been no pharmacogenomic studies for T2D in which a population of color was used in the discovery analysis. This study sought to identify genomic variants associated with metformin response in African American patients with diabetes. RESEARCH DESIGN AND METHODS: Patients in the discovery set were adult, African American participants from the Diabetes Multi-omic Investigation of Drug Response (DIAMOND), a cohort study of patients with T2D from a health system serving southeast Michigan. DIAMOND participants had genome-wide genotype data and longitudinal electronic records of laboratory results and medication fills. The genome-wide discovery analysis identified polymorphisms correlated to changes in glycated hemoglobin (HbA1c) levels among individuals on metformin monotherapy. Lead associations were assessed for replication in an independent cohort of African American participants from Kaiser Permanente Northern California (KPNC) and in European American participants from DIAMOND. RESULTS: The discovery set consisted of 447 African American participants, whereas the replication sets included 353 African American KPNC participants and 466 European American DIAMOND participants. The primary analysis identified a variant, rs143276236, in the gene ARFGEF3, which met the threshold for genome-wide significance, replicated in KPNC African Americans, and was still significant in the meta-analysis (P = 1.17 × 10-9). None of the significant discovery variants replicated in European Americans DIAMOND participants. CONCLUSIONS: We identified a novel and biologically plausible genetic variant associated with a change in HbA1c levels among African American patients on metformin monotherapy. These results highlight the importance of diversity in pharmacogenomic studies.

Cardiac contraction and relaxation are regulated by distinct subcellular cAMP pools.

Cells interpret a variety of signals through G-protein-coupled receptors (GPCRs) and stimulate the generation of second messengers such as cyclic adenosine monophosphate (cAMP). A long-standing puzzle is deciphering how GPCRs elicit different physiological responses despite generating similar levels of cAMP. We previously showed that some GPCRs generate cAMP from both the plasma membrane and the Golgi apparatus. Here we demonstrate that cardiomyocytes distinguish between subcellular cAMP inputs to elicit different physiological outputs. We show that generating cAMP from the Golgi leads to the regulation of a specific protein kinase A (PKA) target that increases the rate of cardiomyocyte relaxation. In contrast, cAMP generation from the plasma membrane activates a different PKA target that increases contractile force. We further validated the physiological consequences of these observations in intact zebrafish and mice. Thus, we demonstrate that the same GPCR acting through the same second messenger regulates cardiac contraction and relaxation dependent on its subcellular location.

SLCO1B1 functional variants and statin-induced myopathy in people with recent genealogical ancestors from Africa: a population-based real-world study.

Background: Clinical pharmacogenetic implementation guidelines for statin therapy are derived from evidence of primarily Eurocentric study populations. Functional SLCO1B1 variants that are rare in these study populations have not been investigated as a determinant of statin myotoxicity and are thus missing from guideline inclusion. Objective: Determine the relationship between candidate functional SLCO1B1 variants and statin-induced myopathy in people with recent genealogical ancestors from Africa. Design: Population-based pharmacogenetic study using real-world evidence from electronic health record-linked biobanks. Setting: Various health care settings. Participants: Self-identified white and Black statin users with genome-wide genotyping data available. Measurements: Primarily, the odds of statin-induced myopathy + rhabdomyolysis. Secondarily, total bilirubin levels. Thirdly, cell-based functional assay results. Results: Meta-analyses results demonstrated an increased risk of statin-induced myopathy + rhabdomyolysis with c.481+1G>T (odds ratio [OR] = 3.27, 95% confidence interval [CI] 1.43-7.46, P =.005) and c.1463G>C (OR = 2.45, 95% CI 1.04-5.78, P =.04) for Black participants. For White participants, c.521T>C was also significantly associated with increased risk of statin-induced myopathy + rhabdomyolysis (OR = 1.41, 95% CI 1.20-1.67, P =5.4x10 -5 ). This effect size for c.521T>C was similar in the Black participants, but did not meet the level of statistical significance (OR = 1.47, 95% CI 0.58-3.73, P =0.41). Supporting evidence using total bilirubin as an endogenous biomarker of SLCO1B1 function as well as from cell-based functional studies corroborated these findings. Limitations: Data limited to severe statin myotoxicity events. Conclusion: Our findings implicate Afrocentric SLCO1B1 variants on preemptive pharmacogenetic testing panels, which could have an instant impact on reducing the risk of statin-associated myotoxicity in historically excluded groups. Primary Funding Source: National Institutes of Health, Office of the Director - All of Us (OD-AoURP).

Structural and molecular basis of choline uptake into the brain by FLVCR2.

Choline is an essential nutrient that the human body needs in vast quantities for cell membrane synthesis, epigenetic modification, and neurotransmission. The brain has a particularly high demand for choline, but how it enters the brain has eluded the field for over fifty years. The MFS transporter FLVCR1 was recently determined to be a choline transporter, and while this protein is not highly expressed at the blood-brain barrier (BBB), its relative FLVCR2 is. Previous studies have shown that mutations in human Flvcr2 cause cerebral vascular abnormalities, hydrocephalus, and embryonic lethality, but the physiological role of FLVCR2 is unknown. Here, we demonstrate both in vivo and in vitro that FLVCR2 is a BBB choline transporter and is responsible for the majority of choline uptake into the brain. We also determine the structures of choline-bound FLVCR2 in the inward- and outward-facing states using cryo-electron microscopy to 2.49 and 2.77 Å resolution, respectively. These results reveal how the brain obtains choline and provide molecular-level insights into how FLVCR2 binds choline in an aromatic cage and mediates its uptake. Our work could provide a novel framework for the targeted delivery of neurotherapeutics into the brain.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates.

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17ß-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

TICBase: Integrated Resource for Data on Drug and Environmental Chemical Interactions with Mammalian Drug Transporters.

Environmental health science seeks to predict how environmental toxins, chemical toxicants, and prescription drugs accumulate and interact within the body. Xenobiotic transporters of the ATP-binding cassette (ABC) and solute carrier (SLC) superfamilies are major determinants of the uptake and disposition of xenobiotics across the kingdoms of life. The goal of this study was to integrate drug and environmental chemical interactions of mammalian ABC and SLC proteins in a centralized, integrative database. We built upon an existing publicly accessible platform-the "TransPortal"-which was updated with novel data and searchable features on transporter-interfering chemicals from manually curated literature data. The integrated resource TransPortal-TICBase (https://transportal.compbio.ucsf.edu) now contains information on 46 different mammalian xenobiotic transporters of the ABC- and SLC-type superfamilies, including 13 newly added rodent and 2 additional human drug transporters, 126 clinical drug-drug interactions, and a more than quadrupled expansion of the initial in vitro chemical interaction data from 1,402 to 6,296 total interactions. Based on our updated database, environmental interference with major human and rodent drug transporters occurs across the ABC- and SLC-type superfamilies, with kinetics indicating that some chemicals, such as the ionic liquid 1-hexylpyridinium chloride and the antiseptic chlorhexidine, can act as strong inhibitors with potencies similar or even higher than pharmacological model inhibitors. The new integrated web portal serves as a central repository of current and emerging data for interactions of prescription drugs and environmental chemicals with human drug transporters. This archive has important implications for predicting adverse drug-drug and drug-environmental chemical interactions and can serve as a reference website for the broader scientific community of clinicians and researchers.

The Impacts of Slc19a3 Deletion and Intestinal SLC19A3 Insertion on Thiamine Distribution and Brain Metabolism in the Mouse.

The Thiamine Transporter 2 (THTR2) encoded by SLC19A3 plays an ill-defined role in the maintenance of tissue thiamine, thiamine monophosphate, and thiamine diphosphate (TDP) levels. To evaluate the impact of THTR2 on tissue thiamine status and metabolism, we expressed the human SLC19A3 transgene in the intestine of total body Slc19a3 knockout (KO) mice. Male and female wildtype (WT) and transgenic (TG) mice were fed either 17 mg/kg (1×) or 85 mg/kg (5×) thiamine hydrochloride diet, while KOs were only fed the 5× diet. Thiamine vitamers in plasma, red blood cells, duodenum, brain, liver, kidney, heart, and adipose tissue were measured. Untargeted metabolomics were performed on the brain tissues of groups with equivalent plasma thiamine. KO mice had ~two- and ~three-fold lower plasma and brain thiamine levels than WT on the 5× diet. Circulating vitamers were sensitive to diet and equivalent in TG and WT mice. However, TG had 60% lower thiamine but normal brain TDP levels regardless of diet, with subtle differences in the heart and liver. The loss of THTR2 reduced levels of nucleic acid and amino acid derivatives in the brain. Therefore, mutation or inhibition of THTR2 may alter the brain metabolome and reduce the thiamine reservoir for TDP biosynthesis.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates.

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17ß-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

The full spectrum of OCT1 (SLC22A1) mutations bridges transporter biophysics to drug pharmacogenomics.

Membrane transporters play a fundamental role in the tissue distribution of endogenous compounds and xenobiotics and are major determinants of efficacy and side effects profiles. Polymorphisms within these drug transporters result in inter-individual variation in drug response, with some patients not responding to the recommended dosage of drug whereas others experience catastrophic side effects. For example, variants within the major hepatic Human organic cation transporter OCT1 (SLC22A1) can change endogenous organic cations and many prescription drug levels. To understand how variants mechanistically impact drug uptake, we systematically study how all known and possible single missense and single amino acid deletion variants impact expression and substrate uptake of OCT1. We find that human variants primarily disrupt function via folding rather than substrate uptake. Our study revealed that the major determinants of folding reside in the first 300 amino acids, including the first 6 transmembrane domains and the extracellular domain (ECD) with a stabilizing and highly conserved stabilizing helical motif making key interactions between the ECD and transmembrane domains. Using the functional data combined with computational approaches, we determine and validate a structure-function model of OCT1s conformational ensemble without experimental structures. Using this model and molecular dynamic simulations of key mutants, we determine biophysical mechanisms for how specific human variants alter transport phenotypes. We identify differences in frequencies of reduced function alleles across populations with East Asians vs European populations having the lowest and highest frequency of reduced function variants, respectively. Mining human population databases reveals that reduced function alleles of OCT1 identified in this study associate significantly with high LDL cholesterol levels. Our general approach broadly applied could transform the landscape of precision medicine by producing a mechanistic basis for understanding the effects of human mutations on disease and drug response.

The Effect of Trimethoprim on Thiamine Absorption: A Transporter-Mediated Drug-Nutrient Interaction.

Trimethoprim is predicted to inhibit several thiamine transporters, including the primary thiamine intestinal absorptive transporter, ThTR-2, and the hepatic and renal organic cation transporters, OCT1, OCT2, and MATEs. To investigate the effect of trimethoprim on thiamine absorption, studies were conducted in cells, mice, and healthy volunteers and supported by use of real-world data. In a randomized, crossover clinical study, seven healthy volunteers were given a single oral dose of thiamine or thiamine plus trimethoprim, followed by blood sampling. The thiamine area under the curve (AUC) increased with trimethoprim co-administration (P value = 0.031). Similar results were seen in mice. Trimethoprim appeared to act on thiamine absorption through inhibition of hepatic OCT1 as evidenced from its ability to modulate levels of isobutyrylcarnitine and propionylcarnitine, OCT1 biomarkers identified from metabolomic analyses. Real-world data further supported this finding, showing an association between trimethoprim use and higher levels of triglycerides, LDL cholesterol, and total cholesterol, consistent with OCT1 inhibition (P values: 2.2 × 10-16 , 5.75 × 10-7 , and 5.82 × 10-7 , respectively). These findings suggest that trimethoprim increases plasma levels of thiamine by inhibiting hepatic OCT1. Trimethoprim reduced urinary excretion and clearance of biomarkers for OCT2 and MATEs, consistent with inhibition of renal organic cation transporters. This inhibition did not appear to play a role in the observed increases in thiamine levels. This study highlights the potential for drug-nutrient interactions involving transporters, in addition to transporters' established role in drug-drug interactions.

Characterization of P-glycoprotein orthologs from human, sheep, pig, dog, and cat.

The ATP-binding cassette transporter P-glycoprotein (P-gp) limits the oral bioavailability of many drugs. Although P-gp has been well studied in humans and mice, little is known about the substrate specificities of many of its species orthologs. To address this, we performed in vitro analysis of P-gp transporter function using HEK293 cells stably expressing human, ovine, porcine, canine, and feline P-gp. We also employed a human physiologically based pharmacokinetic (PBPK) model to assess variations in digoxin exposure resulting from altered P-gp function. Compared to human P-gp, sheep P-gp had significantly less digoxin efflux (2.3-fold ±0.04 vs. 1.8-fold ±0.03, p < .0001) and all species orthologs had significantly less quinidine efflux compared with human P-gp (p < .05). Human P-gp also had significantly greater efflux of talinolol compared to sheep and dog P-gp (1.9-fold ±0.04 vs. 1.6-fold ±0.06, p = .003 and 1.6-fold ±0.05, p = .0002, respectively). P-gp expression protected all lines against paclitaxel-induced toxicity, with sheep P-gp being significantly less protective. The inhibitor verapamil demonstrated dose-dependent inhibition of all P-gp orthologs. Finally, a PBPK model showed digoxin exposure was sensitive to altered P-gp activity. Overall, our study found that species differences in this major drug transporter exist and that the appropriate species ortholog of P-gp should be evaluated during veterinary drug development.

Identification of Genetic Variation Influencing Metformin Response in a Multiancestry Genome-Wide Association Study in the Diabetes Prevention Program (DPP).

Genome-wide significant loci for metformin response in type 2 diabetes reported elsewhere have not been replicated in the Diabetes Prevention Program (DPP). To assess pharmacogenetic interactions in prediabetes, we conducted a genome-wide association study (GWAS) in the DPP. Cox proportional hazards models tested associations with diabetes incidence in the metformin (MET; n = 876) and placebo (PBO; n = 887) arms. Multiple linear regression assessed association with 1-year change in metformin-related quantitative traits, adjusted for baseline trait, age, sex, and 10 ancestry principal components. We tested for gene-by-treatment interaction. No significant associations emerged for diabetes incidence. We identified four genome-wide significant variants after correcting for correlated traits (P < 9 × 10-9). In the MET arm, rs144322333 near ENOSF1 (minor allele frequency [MAF]AFR = 0.07; MAFEUR = 0.002) was associated with an increase in percentage of glycated hemoglobin (per minor allele, ß = 0.39 [95% CI 0.28, 0.50]; P = 2.8 × 10-12). rs145591055 near OMSR (MAF = 0.10 in American Indians) was associated with weight loss (kilograms) (per G allele, ß = -7.55 [95% CI -9.88, -5.22]; P = 3.2 × 10-10) in the MET arm. Neither variant was significant in PBO; gene-by-treatment interaction was significant for both variants [P(G×T) < 1.0 × 10-4]. Replication in individuals with diabetes did not yield significant findings. A GWAS for metformin response in prediabetes revealed novel ethnic-specific associations that require further investigation but may have implications for tailored therapy.

Initial Insights into the Genetic Variation Associated with Metformin Treatment Failure in Youth with Type 2 Diabetes.

Metformin is the first-line treatment for type 2 diabetes (T2D) in youth but with limited sustained glycemic response. To identify common variants associated with metformin response, we used a genome-wide approach in 506 youth from the Treatment Options for Type 2 Diabetes in Adolescents and Youth (TODAY) study and examined the relationship between T2D partitioned polygenic scores (pPS), glycemic traits, and metformin response in these youth. Several variants met a suggestive threshold (P < 1 × 10-6), though none including published adult variants reached genome-wide significance. We pursued replication of top nine variants in three cohorts, and rs76195229 in ATRNL1 was associated with worse metformin response in the Metformin Genetics Consortium (n = 7,812), though statistically not being significant after Bonferroni correction (P = 0.06). A higher ß-cell pPS was associated with a lower insulinogenic index (P = 0.02) and C-peptide (P = 0.047) at baseline and higher pPS related to two insulin resistance processes were associated with increased C-peptide at baseline (P = 0.04,0.02). Although pPS were not associated with changes in glycemic traits or metformin response, our results indicate a trend in the association of the ß-cell pPS with reduced ß-cell function over time. Our data show initial evidence for genetic variation associated with metformin response in youth with T2D.

The Clinical Pharmacogenetics Implementation Consortium Guideline for SLCO1B1, ABCG2, and CYP2C9 genotypes and Statin-Associated Musculoskeletal Symptoms.

Statins reduce cholesterol, prevent cardiovascular disease, and are among the most commonly prescribed medications in the world. Statin-associated musculoskeletal symptoms (SAMS) impact statin adherence and ultimately can impede the long-term effectiveness of statin therapy. There are several identified pharmacogenetic variants that impact statin disposition and adverse events during statin therapy. SLCO1B1 encodes a transporter (SLCO1B1; alternative names include OATP1B1 or OATP-C) that facilitates the hepatic uptake of all statins. ABCG2 encodes an efflux transporter (BCRP) that modulates the absorption and disposition of rosuvastatin. CYP2C9 encodes a phase I drug metabolizing enzyme responsible for the oxidation of some statins. Genetic variation in each of these genes alters systemic exposure to statins (i.e., simvastatin, rosuvastatin, pravastatin, pitavastatin, atorvastatin, fluvastatin, lovastatin), which can increase the risk for SAMS. We summarize the literature supporting these associations and provide therapeutic recommendations for statins based on SLCO1B1, ABCG2, and CYP2C9 genotype with the goal of improving the overall safety, adherence, and effectiveness of statin therapy. This document replaces the 2012 and 2014 Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines for SLCO1B1 and simvastatin-induced myopathy.

High Throughput Screening of a Prescription Drug Library for Inhibitors of Organic Cation Transporter 3, OCT3.

INTRODUCTION: The organic cation transporter 3 (OCT3, SLC22A3) is ubiquitously expressed and interacts with a wide array of compounds including endogenous molecules, environmental toxins and prescription drugs. Understudied as a determinant of pharmacokinetics and pharmacodynamics, OCT3 has the potential to be a major determinant of drug absorption and disposition and to be a target for drug-drug interactions (DDIs). GOAL: The goal of the current study was to identify prescription drug inhibitors of OCT3. METHODS: We screened a compound library consisting of 2556 prescription drugs, bioactive molecules, and natural products using a high throughput assay in HEK-293 cells stably expressing OCT3. RESULTS: We identified 210 compounds that at 20 µM inhibit 50% or more of OCT3-mediated uptake of 4-Di-1-ASP (2 µM). Of these, nine were predicted to inhibit the transporter at clinically relevant unbound plasma concentrations. A Structure-Activity Relationship (SAR) model included molecular descriptors that could discriminate between inhibitors and non-inhibitors of OCT3 and was used to identify additional OCT3 inhibitors. Proteomics of human brain microvessels (BMVs) indicated that OCT3 is the highest expressed OCT in the human blood-brain barrier (BBB). CONCLUSIONS: This study represents the largest screen to identify prescription drug inhibitors of OCT3. Several are sufficiently potent to inhibit the transporter at therapeutic unbound plasma levels, potentially leading to DDIs or off-target pharmacologic effects.

Emerging Roles of the Human Solute Carrier 22 Family.

The human Solute Carrier 22 family (SLC22), also termed the organic ion transporter family, consists of 28 distinct multi-membrane spanning proteins, which phylogenetically cluster together according to their charge specificity for organic cations (OCTs), organic anions (OATs) and organic zwitterion/cations (OCTNs). Some SLC22 family members are well characterized in terms of their substrates, transport mechanisms and expression patterns, as well as their roles in human physiology and pharmacology, whereas others remain orphans with no known ligands. Pharmacologically, SLC22 family members play major roles as determinants of the absorption and disposition of many prescription drugs, and several including the renal transporters, OCT2, OAT1 and OAT3 are targets for many clinically important drug-drug interactions. In addition, mutations in some of these transporters (SLC22A5 (OCTN2) and SLC22A12 (URAT1) lead to rare monogenic disorders. Genetic polymorphisms in SLC22 transporters have been associated with common human disease, drug response and various phenotypic traits. Three members in this family were deorphaned in very recently: SLC22A14, SLC22A15 and SLC22A24, and found to transport specific compounds such as riboflavin (SLC22A14), anti-oxidant zwitterions (SLC22A15) and steroid conjugates (SLC22A24). Their physiologic and pharmacological roles need further investigation. This review aims to summarize the substrates, expression patterns and transporter mechanisms of individual SLC22 family members and their roles in human disease and drug disposition and response. Gaps in our understanding of SLC22 family members are described. Significance Statement In recent years, three members of the SLC22 family of transporters have been deorphaned and found to play important roles in the transport of diverse solutes. New research has furthered our understanding of the mechanisms, pharmacological roles, and clinical impact of SLC22 transporters. This minireview provides overview of SLC22 family members of their physiologic and pharmacologic roles, the impact of genetic variants in the SLC22 family on disease and drug response, and summary of recent studies deorphaning SLC22 family members.

Genome-Wide Meta-analysis Identifies Genetic Variants Associated With Glycemic Response to Sulfonylureas.

OBJECTIVE: Sulfonylureas, the first available drugs for the management of type 2 diabetes, remain widely prescribed today. However, there exists significant variability in glycemic response to treatment. We aimed to establish heritability of sulfonylurea response and identify genetic variants and interacting treatments associated with HbA1c reduction. RESEARCH DESIGN AND METHODS: As an initiative of the Metformin Genetics Plus Consortium (MetGen Plus) and the DIabetes REsearCh on patient straTification (DIRECT) consortium, 5,485 White Europeans with type 2 diabetes treated with sulfonylureas were recruited from six referral centers in Europe and North America. We first estimated heritability using the generalized restricted maximum likelihood approach and then undertook genome-wide association studies of glycemic response to sulfonylureas measured as HbA1c reduction after 12 months of therapy followed by meta-analysis. These results were supported by acute glipizide challenge in humans who were naïve to type 2 diabetes medications, cis expression quantitative trait loci (eQTL), and functional validation in cellular models. Finally, we examined for possible drug-drug-gene interactions. RESULTS: After establishing that sulfonylurea response is heritable (mean ± SEM 37 ± 11%), we identified two independent loci near the GXYLT1 and SLCO1B1 genes associated with HbA1c reduction at a genome-wide scale (P < 5 × 10-8). The C allele at rs1234032, near GXYLT1, was associated with 0.14% (1.5 mmol/mol), P = 2.39 × 10-8), lower reduction in HbA1c. Similarly, the C allele was associated with higher glucose trough levels (ß = 1.61, P = 0.005) in healthy volunteers in the SUGAR-MGH given glipizide (N = 857). In 3,029 human whole blood samples, the C allele is a cis eQTL for increased expression of GXYLT1 (ß = 0.21, P = 2.04 × 10-58). The C allele of rs10770791, in an intronic region of SLCO1B1, was associated with 0.11% (1.2 mmol/mol) greater reduction in HbA1c (P = 4.80 × 10-8). In 1,183 human liver samples, the C allele at rs10770791 is a cis eQTL for reduced SLCO1B1 expression (P = 1.61 × 10-7), which, together with functional studies in cells expressing SLCO1B1, supports a key role for hepatic SLCO1B1 (encoding OATP1B1) in regulation of sulfonylurea transport. Further, a significant interaction between statin use and SLCO1B1 genotype was observed (P = 0.001). In statin nonusers, C allele homozygotes at rs10770791 had a large absolute reduction in HbA1c (0.48 ± 0.12% [5.2 ± 1.26 mmol/mol]), equivalent to that associated with initiation of a dipeptidyl peptidase 4 inhibitor. CONCLUSIONS: We have identified clinically important genetic effects at genome-wide levels of significance, and important drug-drug-gene interactions, which include commonly prescribed statins. With increasing availability of genetic data embedded in clinical records these findings will be important in prescribing glucose-lowering drugs.

Genomewide Association Studies in Pharmacogenomics.

The increasing availability of genotype data linked with information about drug-response phenotypes has enabled genomewide association studies (GWAS) that uncover genetic determinants of drug response. GWAS have discovered associations between genetic variants and both drug efficacy and adverse drug reactions. Despite these successes, the design of GWAS in pharmacogenomics (PGx) faces unique challenges. In this review, we analyze the last decade of GWAS in PGx. We review trends in publications over time, including the drugs and drug classes studied and the clinical phenotypes used. Several data sharing consortia have contributed substantially to the PGx GWAS literature. We anticipate increased focus on biobanks and highlight phenotypes that would best enable future PGx discoveries.